The chiral condensate in a constant electromagnetic field at 0{p^) 
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We examine tlie shift in the chiral condensate due to a constant electromagnetic field at C(p^) 
using SU{2) chiral perturbation theory and a realistic — 140 MeV. We find that this value 
differs significantly fi'om the value calculated using Mtt — 0, while the magnitude of the two-loop 
correction is unclear due to the uncertainty in the experimentally determined value of the relevant 
£6 LEG. 
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I. INTRODUCTION 

In QCD, the chiral condensate is important as it is 
the order parameter of chiral symmetry breaking. As 
such, its behavior is key in understanding QCD in ex- 
treme conditions. Finite temperature and pressure ef- 
fects have been studied extensively. On the other hand, 
comparatively few efforts have focused on the effects of fi- 
nite electromagnetic fields. Prior studies have been done 
either using chiral perturbation theory (xPT)[l|,[3,Q, the 
effective theory for low-energy QCD, or models compati- 
ble with the large- iVc expansion0,Q (the NJL model, in 
particulardHHIil). The method used in this paper is 
XPT. This approach has the advantage of being model- 
independent and systematic, but has the disadvantage of 
containing a number of undetermined parameters (low- 
energy constants or LECs). We extend prior work done 
at the = limit 0, |3| and at one loop 10] to two 
loops for = 140 MeV. 

Other calculations using xPT have focused on A/^ = 0, 
and while this may be of interest theoretically, it has 
at best a narrow window of validity p^. In addition, it 
is unlikely to find in nature a real electric or magnetic 
field with eE » or eH Mj, which is required 
for the approximation to be reasonable. Incidentally, the 
opposite limit, ^ eH, might be of more interest in 
the sense that it can be produced in the laboratory, but 
in such a regime the shift in the condensate is miniscule. 

Continuing a calculation to higher orders in an expan- 
sion is always of interest at least in the trivial sense of 
finding a more precise result. In this case, large-A^c 
QCD[ii|, provides another possible motivation for 
why the 0{p^) result might be of interest. Low-energy 
constants (LECs) of the same chiral order will have dif- 
ferent orders oi Nc, depending upon the number of flavor 
traces in the term they multiply. This can be understood 
as follows: a flavor trace corresponds to a quark loop in 
the analogous QCD calculation, and large-A^c counting 
rules indicate quark loops are down by a power of Nc, 
and the LEC is the only parameter available to absorb 
this difference. 

This large-iV(7 dependence of the LECs can provide 
hints to the convergence of the chiral expansion. There 



are several processes which have been calculated at two- 
loop order, most of which show a close match between 
the 0{p^) results and the experimental results p^ . 
The 0{p^) correction required to achieve this agreement 
varies, but in particular, the 0{p^) calculation of the pro- 
cess 77 — > TT^i" is strikingly more accurate than the 
0(p4) result [m]. This is also the process which provides 
a (very rough) estimate of the LEC we use here, so it is 
conceivable that the 0{p^) correction will be important 
in our case, as well. 

We will first proceed with a brief overview of SU{2) 
xPT. We then follow with an analytical calculation of 
the shift in the chiral condensate due to a magnetic field 
at 0{p^ ), and finally, a numerical analysis of the shift for 
general electromagnetic fields. 



II. CHIRAL PERTURBATION THEORY 
A. Basics 

This section is a brief summary of the chiral pertur- 
bation theory notation that will be used in this paper. 
For a detailed description of the theory, see the origi- 
nal papers by Gasser and LeutwylerQ or a number of 
reviews (refs. [H, [l^, for example). Treatment of the 
chiral lagrangian to 0{p^) can be found in [Tsl. [TtI [T§|. 

The building blocks that are used to construct the chi- 
ral lagrangian include U — u{(f>)'^ , containing the dynam- 
ical pion fields, and external fields s, p, a^, and v^. We 
will consider here SU (2) flavor symmetry with m„ « rrid, 
with external fields corresponding to only a constant elec- 
tromagnetic field and quark masses. We then have 



X = 2B{s + ip) = 2BM 



A., 



(1) 



We also need to define the covariant derivative as 



(2) 
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From these, we define the following operators that will 
contribute to the terms relevant in this paper (in the 



2 



general SU{n) notation): 



(3) 



where are defined 



X'u 



(4) 



Note that in the case we are considering, F^"^ — F'^^ . 

U can be paramaterized in several ways, but we will 
use the Weinberg parameterization 



U = a 



F2 



= 1. 



(5) 



Here, tt" are still the dynamical fields, and a is repre- 
sented as an expansion in terms of tt" from the second 
equation. 

Using these definitions, £2 and £4 can be written as 
follows (in the general SU{N) form), where (A) denotes 
the trace of A 1141: 



£2 ~{u^u^ + x+) 



il. 



(6) 
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contact terms 



The calculation here is up to O(p^), which means that 
we will be using the £2 lagrangian up to two loops and 
the £4 lagrangian up to one loop; the lagrangian will 
also contribute at tree level. This lagrangian has been 
calculated in ref. [T8| . and has more terms than we will 
list (112 for SU(n) and 53 for S'C/(2)). Fortunately, only 
one of these (in SU{2)) will be relevant for our calcula- 
tion, as we will see later, and it can be expressed as: 



(7) 



B. Renormalization 



Renormalization of the theory to 0{p^) was calculated 
in ref. [l|. It has also more recently been calculated for 
the 0{p^) lagrangian in p7j . 

Using the SU (2) LECs, the renormalized couplings can 



be written: 



1 

167r2(d-4) 
1 2 

71 = g,72 g,73 



1 



:,74 = 2, 



(8) 



1 

75 = ~e'76 



1 



r,77 = 



When M ^ 0, these can be expressed in terms of scale- 
independent parameters as: 



(9) 



The renormalization of the Cq term that we will be using 
later can be expressed similarly in terms of the renormal- 
ized LECs from £4 as 



Ci 



(c//) 



2((i-4) 



F^ 



ir , ^ir (1) (2) n 

734 = -^5 + 2^6: 734 = 734 = 



(cKm, d) - 7f ^A^ - (7f ) + 7^ V, d))A) 

(10) 



III. CALCULATION OF E FROM VACUUM 
ENERGY 

The term in the QCD lagrangian which is relevant to 
E is rriqqq. Since we know that S (qq), in the isospin 
limit of mu = nid = m we can calculate the condensate 
from the vacuum energy as 



de 



dm 



(11) 



To first order in M^/F^, the Gell-Mann-Oakes-Renner 
relation F^M^ = S](m„ -I- TOd)[i^ applies, and can be 
used to calculate the shift in the condensate. Unfortu- 
nately, we will need the next order result, which will in- 
troduce some ambiguity as follows. Neither the chiral 
condensate or the quark mass can be defined indepen- 
dently; only their product has a physical meaning. In 
the language of xPT[i|, 



2mE = F^M^ i 1 



Ml 



32^2^2 



(4/ii - [3) + 0{Mt) 



(12) 

The ambiguity here is codified in the unphysical LEG /ii, 
which will vary according to the renormalization conven- 
tion. This is an ambiguity in the definition of E. In order 
to avoid this difficulty, we choose to normalize our results 
according to the quantity Eg, which we define by 



2toEo 



(13) 



with F^ and at their physical values. We will thus 
express our results in terms of AE/Eq, where AE = 



3 




4) 'VWW\^^5)'WWW 



(c) 



(d) 



(e) 



FIG. 1: Diagrams contributing to the vacuum energy shift 
due to an electromagnetic field. Dashed lines denote tt" and 
solid lines denote tt^. 



— I]{H = 0). As we will see, this ratio is unam- 
biguous at the order to which we work. 

We will also need the relationships between M and 
as well as F and (M^ = 2Brh). The differ- 
ence between the lowest-order 0{p'^) result and the 0{p^) 
result is Oip^); we therefore only need one order of 
corrections [l|: 



1 



M2 



+ 0{M^ 



(14) 



rud up to the relevant terms in £g 's'l , including only terms 
which will contribute to the diagrams in Figure [TJ 



2^ M 7 2 



1 



1 2 



(15) 



£4 = 



£6 =4c34Af'(eF^,)2 

Here, terms proportional to and I4 have been added to 
the lagrangian from ref . Q , which are down by an order 
of but have the same overall chiral order. 

As a first case, we will work with the case of pure mag- 
netic fields, where {eF^^)^ = 2{eH)^. This simplifies the 
calculations and allows us to obtain an analytic result. 
We will later generalize to an arbitrary combination of E 
and H fields for numerical analysis. 

The propagator for a scalar particle in a constant H 
field was first calculated in ref. [231, and here we use the 
convenient form also used in ref. |3|: 



This will only be applicable in the first term of the ex- 
pansion, which is two powers of momentum less than the 
maximum order for the calculation. 

Our object is to calculate the condensate for the case of 
a constant electromagnetic field at 0{p^) in xPT. As de- 
scribed above, we can accomplish this by calculating the 
vacuum energy to the same order. As noted by ref. [11, 
this calculation will involve two-loop diagrams with C2 
vertices, one-loop diagrams with an £4 vertex, and tree- 
level diagrams with an £g vertex. Diagrams which con- 
tribute to the vacuum energy will contain only external 
photon lines coming from the constant EM field (which 
is the "vacuum" in this case). Insertions of the electro- 
magnetic field at the £2 level are calculated as part of the 
propagator of the tt*. Thus, the 0{p*) calculation cor- 
responds roughly to a closed single propagator and was 
calculated in ref. 10]. Contributing diagrams must be 
dependent on the electromagnetic field, either through a 
direct insertion or through the propagator of the tt^ . 

With these criteria, we find that the diagrams con- 
tributing to our calculation are as pictured in Fig.[TJ Re- 
turning, then, to the xPT lagrangian, we find that only 
the terms proportional to I3, I4, I5, and Iq can contribute 
from the £4 lagrangian, and only the term proportional 
to C34 (as we anticipated above) can contribute to the 
vacuum energy from Cq. 

We simplify the chiral lagrangians for SU{2) and rn„ — 



D"{k) 



Jk{x-y) p,H 



ds 



cosh{eHs) 



(16) 



where $ 



kl 



kl 



{x,y) = exp{iejy A^{z)dz^,}, kf^ ^ kj + 



kl and 



We will also need the scalar propagator 

d'^k 



D{0) =D{x,x) = J 



A 



=2M2(cAt) 



167r2(d-4)' 



1 



A 



32n' 



(17) 



D{0) and (0) = D^{x,x) are both divergent quanti- 
ties, whereas (0) — D{0) is finite: 



D^"{0) =D"{0)-D{0) 
eH 

X 



167r2 



00 7 

( 1 - 







sinhx 



(18) 



with P — M /eH. This is the same integral as was cal- 
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culated in ref. and can be expressed analytically as: 
eH 



D^" (0) 



/H(/3)=log(2^) + /31og( ^ 



^- 2iogr 



1 + /3 



(19) 



With these, the diagrams in Fig. [T] can be calculated 
fairly straightforwardly to be 0] 



(2) 



A/' 



1(a) 2F2 



i2) 



1(b) 



(2) 



1(C) 



-(2) 



F2 

M4 



(2^5 



2{eHf 



(20) 



1(d) 



1(e) 



=2^3^15^(0) 



The only one of these diagrams which is not expressed 
solely in terms of (0) and -D(O) is q^^^ [, which van- 



ishes generally as well as in the A/jr = case. 

In order to make the divergences and scale-dependence 
explicit, we first make the substitution D^{0) = 
D'^^ {0) + D{0). Any term which is dependent on neither 
H nor D^{0) can then be re-absorbed into the vacuum 
energy; we are only looking for the shift due to H. After 



J 



(2) 

making this substitution, we see that q\ I ^ | is divergent 

11(a) 



and cancelled by a counterterm generated by the com- 

,(2) 



bination (2/3 -|- ^^J^ has both a finite piece, which 
will contribute to the calculation, and a divergent piece. 



-(2) 



(2) 

which is cancehed by a counterterm in C34. / ^^ | and 



1(d) 



are finite, aside from the aforementioned countert- 



1(e) 
eims 

Combining, then, all of these terms, we find the vac- 
uum energy to be 



.(2) 



{eH) 



(167r2)2F2 



l3d{M^ 



(21) 



where we have defined the scale-independent quantity 

1 



M2 

/i2 



(22) 



We then substitute Eq. ^ to find from AP in the 
first -order term. We find that the I3 term cancels, and 
that the correction to F does not play a role (as it only 
appears at the second order). Then, taking a deriva- 
tive, and applying the Gell-Mann-Oakes-Renner relation 
as above, with the first-order corrections for the ©(p^) 
term, we find (/3^ = M^/eH): 



AE(iJ) eH 



Ih{M + 



eH 



167r2F2 



-3(^6-^5) 



r 



1 + log 2 + V 



d{eH) 



(23) 



with ip{x) = ^logr(x) 

Taking — > 0, ip {^) = — 7e, and the shift we find 
agrees with the expression found in ref. [Sj] for the case 
M„ = 0. 

For the case of an E field, we can make the substitution 
H iE and get a similar analytic expression. The E ■ 
H ^ case is somewhat more complicated, but we will 
write an integral expression which wc can later evaluate 
numerically. 

For a convenient paramaterization of the general case, 



J 



we introduce the variables (f) and / such that with J- 
Sl-Ml^lF^^^ndg^E-H^, 



/2 cos(20) 



f sin(20) 



(24) 



Expressed in terms of these variables, the shift in the 
condensate due to an arbitrary combination of fields will 
become (/?/ = Ml/ef): 



AI](.F, g) 
So 

lEHiPfA) 



ef 



167r2F2 
dz 



1 - 



167r2F2 

z2 sin 20 



cos20<{ -([g- [5) (/^^ (/?/,, 



l) + d{Ml] 



2 sin(2; sin (j)) sinh(z cos ( 



le 



(25) 
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This is the same integral as in ref. and as be- 

fore we have had to avoid some potential ambiguity. The 
poles in the integrand indicate an instability in the sys- 
tem, which is interpreted as due to pair creation in an 
electric field We have chosen to regulate the diver- 
gence in a manner which has an imaginary part corre- 
sponding to this pair creation. The magnitude of the 
imaginary part indicates the importance of this instabil- 
ity, though some caution is warranted in interpreting it 
quantitatively. This issue was discussed in more detail in 
ref. 53]. 

Equations (|23)) and (p5)) are the principal results of this 
work. 



IV. NUMERICAL RESULTS 

Because the behavior of the theory is encoded in the 
LECs, a real-world interpretation of the low-energy be- 
havior requires the use of measured LECs. For the £4 
LECs, this is straightforward, as these are individually 
determined with relatively small error. On the other 
hand, £g LECs, such as C34, are more problematic, as 
there are in general many more LECs than easily measur- 
able processes to determine them. These LECs are often 
estimated (at a particular scale) by resonance exchange. 
Unfortunately, C34 in particular is difficult to extract, as 
the resonance processes to which it contributes involve 
only scalar exchange, and because it appears squared in 
these processes, its sign is undetermined. This resonance 



exchange occurs at a scale Mp = 768 [IJ], and it is the 



(scale-dependent) value determined by experiment that 
has an undetermined sign. The scale-independent d is 
positive in both cases. 

The values we use for these constants are VA uM : 




FIG. 2: A comparison of the shift due to a pure magnetic 
field in the KU = case to the A/^ = 140 MeV case. Shaded 
regions indicate uncertainty due to the Cq constant (T . 



Using this method, we plot the total shift in the con- 
densate up to 0{p^) from a general E and H field in 
Fig. [3l and in Fig. [H we plot the ratio of the added cor- 
rection at 0{p^) to the total shift. In these plots, we have 
included a shaded region to indicate the possible values 
for the shift based on a range for (F of (—3,3). 

The asymptotic expression for the shift as (3^^ ^ 00 (for 
an H field, which also provides some qualitative insight 
to other cases) is 



AS(£) 

So 



eH ( Fl 



h — h 
487r2 



167r2 



(28) 



;6-^5 = 3.0±0.3 
d'^ = 8(167r2)2c54 = ±L5±1.5. 



(26) 



With these experimental values, we can plot realistic 
values of the shift in the condensate. In Fig. [2l for the 
case of a pure magnetic field, we compare the value calcu- 
lated for a finite to that for = 0. It is clear that 
these values are significantly different, as in the 0{p'^) 
casejT^]. In this and the figures following, we have chosen 
to extend our results up to ef ~ 290 MeV (the expansion 
parameter is A = AttFt^ = 1.2 GeV). 

Here, we use the same numerical trick as in ref. (To| 
to extract the principal value of the integral numerically. 
We remove the singularities due to the poles located at 
Zi with residue Ri by subtracting the expression 



1 



1 



Z + Zn 



(27) 



The principal value of the integral of this expression is 
zero, but it has a singularity at Zi which exactly cancels 
the singular behavior of the integrand in Ieh- 



which is, of course, zero for 00 {H — *■ 0). This 

expression encodes low-energy behavior for a more real- 
istic regime,; namely, that of the actual pion mass and a 
very small magnetic field. We see from Fig. [4] and in ([28]) 
that the (unknown) sign of has a profound impact on 
the importance of the 0{p^) calculation. The shift in the 
condensate for a positive d^ is significant enough even as 
/ ^ 0, whereas the shift for a negative one is negligible 
up to large values of ef /F^. 

Another notable feature is that the contribution to the 
imaginary part is larger at 0{p^) order as a pure H field 
is approached (while the total imaginary part is going 
to zero). Except in this case where the total imaginary 
part is negligible, the fraction of the imaginary part at 
two loops will be much less significant than its real coun- 
terpart in regimes where the chiral expansion would be 
expected to converge (ef Ml or below). 

The calculation as a whole will only be valid when 
the real part is significantly larger than the imaginary 
part. When the imaginary part dominates, the system 
will break down due to the instability from pair creation. 
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2 4 6 8 10 F,^ 

(b) 

FIG. 3: The imaginaxy and real parts of the total value of the 
shift in the condensate due to general E and H fields, with 
/ and as defined in the text. Shading depicts uncertainty 
due to C34. 



V. DISCUSSION AND CONCLUSIONS 

We have studied the shift in the chiral condensate due 
to an electromagnetic field using chiral perturbation the- 
ory, which is a powerful tool for analyzing the low-energy 
behavior of QCD. Our analysis was done at 0{p^') with 
= 140 MeV. It is obvious that the inclusion of a 
nonzero pion mass greatly affects the result. 

The importance of the 0{p^) correction is less clear. 
Large- iVc reasoning coupled with the results of model- 
based calculations give circumstantial evidence that it 



could play an important role in the final result. How- 
ever, because the sign of the relevant LEG at Cq is un- 
determined by experiment, its effect at 0{p^) could be 
significant or virtually irrelevant. 




-0.5 - 



-1.0 L 

(a) 




(b) 

FIG. 4: The imaginary and real parts of the ratio of the shift 
at two loops to the total shift for the case of general E and 
H fields, with / and 4> as defined in the text. Shading depicts 
uncertainty due to C34. 
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